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Abstract 

We analyze single particle momentum spectra and interferometry radii in central Au+Au 
collisions at RHIC within a hydro-inspired parametrization accounting for continuous hadron 
emission through the whole lifetime of hydrodynamically expanding fireball. We found that 
a satisfactory description of the data is achieved for a physically reasonable set of parameters 
when the emission from non space-like sectors of the enclosed freeze-out hypersurface is fairly 
long: 9 fm/c. This protracted surface emission is compensated in outward interferometry radii 
by positive rout — t correlations that are the result of an intensive transverse expansion. The 
main features of the experimental data are reproduced: in particular, the obtained ratio of 
the outward to sideward interferometry radii is less than unity and decreases with increasing 
transverse momenta of pion pairs. The extracted value of the temperature of emission from the 
surface of hydrodynamic tube approximately coincides with one found at chemical freeze-out in 
RHIC Au-|-Au collisions. A significant contribution of the surface emission to the spectra and 
to the correlation functions at relatively large transverse momenta should be taken into account 
in advanced hydrodynamic models of ultrarelativistic nucleus-nucleus collisions. 
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1 Introduction 

It is quite possible that a new state of matter is created at unprecedented high energy density 
reached in relativistic heavy-ion collisions at RHIC This conclusion is based, in particular, 
on the success of hydrodynamic models in description of the measured single-particle hadron mo- 
mentum spectra and elliptic flows Within the hydrodynamic approach, there is, in principle, 
a possibility to extract the equation of state (EoS) of the thermalized system and then get infor- 
mation about the early stage of the collision processes and even the type of the phase transition 
between QGP and the hadron matter. 

However quantitative determination of the EoS in the hydrodynamic approach is a nontrivial 
problem. It depends on both the initial conditions and system "freeze-out" prescription. The 
latter determines the final conditions of hydrodynamic expansion: one cannot use hydrodynamic 
equations for infinitely large times since the resulting very small densities destroy the picture of 
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continuous medium of real particles. Usually, the freeze-out hypersurface which confines the 4- 
volume of hydrodynamic evolution is included as an external input with respect to hydrodynamic 
equations. Then, typically, particle spectra are described by the Cooper-Frye prescription (CFp) (3] 
which treats freeze-out as a sudden transition from local thermal equilibrium to free streaming that 
happens on some space-like 3D hypersurface, e.g., on a space-like part of the isotherm: T ~ m^r [1]. 
In reality the process of freeze-out, or particle liberation, is quite complicated because the particles 
escape from the system during the whole period of its evolution (see, e.g., p]). A method which is 
based on Boltzmann equation was recently proposed in Ref. [Sj to describe the spectra formation 
in the hydrodynamic approach to A+A collisions using the escape probabilities. The region of 
applicability of CFp and its possible generalizations are discussed there, as well as in Ref. [7]. In 
present work we do not consider the whole complexity of the freeze-out process (see, e.g., Refs. 
[nilZllHlin])- Instead we account for continuous (in time) emission process by applying the generalized 
Cooper-Frye prescription for an enclosed freeze-out hypersurface. Such an enclosed hypersurface 
contains space-like sectors of the volume emission and non-space-like ones (with space-like normal) 
of the surface emission. 

It is well known ^JJl that protracted surface emission may lead to the ratio of the outward to 
sideward interferometry radii, Rout/Rside-, much bigger than unity for typical scenarios of phase 
transition in heavy ion collisions These expectations are, however, in contradiction with 

current experimental data from RHIC where Rout/Rside ~ 1 (see, e.g., ^H). This discrepancy is 
a component of the so-called "HBT puzzle" JHl- Naively, from this one might conclude that the 
duration of emission is very short. However, the Rout/Rside ratio is also quite sensitive to the shape 
of freeze-out hypersurface because, unlike Rside, Rout is a mixture of the out- width, time spread of 
emission, and the rout — t correlations. Indeed, in the Gaussian approximation for the correlation 
function the interferometry radii can be expressed in terms of the space-time variances^ (see, e.g., 

my- 

Rfip) = ((Ar, - v^Atf)p = (Arf), + vf{At\ - 2v^{AriAt),, (1) 

where Vi = Pi/po, Pi, fi [i = out, side, long) are the Cartesian components of the vectors v, p, 
and r, respectively, and p^ = {po,Pout,0,Piong) is the mean 4-momentum of the two registered 
particles. Here Ar^ = r^ — {ri)p, At = t — {t)p, and {■■.)p denotes the averaged (over the distribution 
function) value taken at some momentum p. Note that in the Bertsch-Pratt reference frame |15j 
Pside = and therefore pout = Pt, where px is the absolute value of the transverse component of 
the vector p. It is easy to see from Eq. (0) that positive rout — t correlations, {AroutAt)p > 0, give 
a negative contribution to Rout reducing thereby the Rout/Rside ratio. ^ Therefore one can conclude 
that relatively small Rout/Rside — 1 ratio in a case of a prolonged surface emission can take place 
if there are positive rout — t correlations in the corresponding sector of the freeze-out hypersurface. 
Such a sector has typically a space-like 4-normal. Then an interval between space-time points in 
such a sector can be, generally, space-like as well as time-like. While the bulk of hadrons are emitted 
from the sector with a time-like normal (volume emission), the surface emission from the sector 
with a space-like normal can be, nevertheless, important for spectra and HBT radii formation'^, 
especially at relatively high px- A difference in temperature values and intensities of collective 

^Note that if the distribution function has a non-Gaussian form, which is typical for expanding sources, then 
the expression of the interferometry radii through space-time variances can be associated with the behavior of the 
correlation function at small q = pi — P2 only, and, therefore, could be used solely for qualitative estimates. 

^This mechanism of reduction of the Rout/Rside ratio is realized in a multiphase transport (AMPT) model 

^See the relevant analysis in Ref. |17| 



2 



flows^ in different freeze-out sectors can increase the effect, and thus a consistent description of 
the single-particle spectra and interferometry radii can be reached in hydrodynamically-motivated 
models with an enclosed freeze-out hypersurface that accounts for continuous (in time) character 
of particle emission. 

It is well known that when the freeze-out hypersurface a has non-space-like sectors, the CFp 
is inconsistent. Then the CFp should be modified to exclude formally negative contributions 
to the particle number at the corresponding momenta. The simplest prescription is to present 
the distribution function as a product of a local thermal distribution and the step function like 
9{p^n^{x)) |19 | l2f) j . n^n^ = if; where is a time- like or space- like outward normal to a freeze-out 
hypersurface a.^ Thereby freeze-out is restricted to those particles for which p^n^^{x) is positive. 

Sometimes, to avoid the above mentioned problems with CFp, the continuity of particles emis- 
sion in heavy ion collisions is taken into account by means of the emission (source) function S{x,p) 
(see, e.g., ^M.), which is used instead of the distribution (Wigner) function f{x,p) in some hydro- 
inspired parameterizations and is usually chosen to be proportional to the local equilibrium dis- 
tribution function, fi.eqix,p), with a smearing (proper) time factor exp(— (r — To)^/Ar^). First, 
such a prescription loses an extremely important information about r — r j correlations: it is nat- 
urally that at early times the emission function is concentrated mostly at the periphery of the 
system on the boundary with vacuum, since particles cannot escape from the hot and dense central 
region. The approximation of surface emission is, thus, much more realistic for the early times. 
Second, while the distribution function is well defined for a given state of the system, like the 
local Bose-Einstein or Fermi-Dirac distributions are associated with a locally equilibrated state, 
S{x,p) is much more complicated and model dependent [Zll^- It is demonstrated in Ref. [S], 
starting from a particular exact solution of the Boltzmann equation, that the emission function is 
not proportional to the local equilibrium distribution (Wigner function), although the system is in 
a local equilibrium state. Therefore the widely used (see, e.g., so called "blast-wave" model j22| ) 
"local equilibrium" ansatz, S ~ fi,eq., is supported neither by analytic kinetic results nor numerical 
transport calculations except for the formal case when all the particles are emitted simultaneously: 
the emission function is then the distribution function multiplied by a delta function concentrated 
on the freeze-out hypersurface. Perhaps, too large values of maximal collective velocities, close to 
the speed of light, which are necessary to fit the spectra in the blast-wave parametrization, and a 
failure in fitting HBT radii (especially Rside, see Refs. j22j) are caused by an improper description 
of the particle emission process. 

The aim of this work is a coherent description of pion, kaon and proton single-particle momen- 
tum spectra as well as the HBT pion radii in central RHIC Au+Au collisions at ^/s^n = 200 GeV 
by means of a simple hydrodynamically-motivated parametrization with an enclosed freeze-out hy- 
persurface. As a matter of fact, we generalize the blast-wave parametrization for a realistic case 
of continuous emission over the whole lifetime of the fireball. We demonstrate that a consistent 
description of the observables is possible if the freeze-out hypersurface has a non-space-like sector 
emitting a noticeable part of the particles.^ The latter condition is satisfied when particles are 
emitted from the surface of an expanding hydrodynamic tube soon after the hadronization, which 

^It was found in Ref. |18| that the transport freeze-out process is similar to evaporation: high-pT particles freeze 
out early from the surface, while 1ow-j3t ones decouple later from the system's center. 

^If a fluid element that crosses a sector of the freeze-out hypersurface with a space-like normal decays, preserving 
its total particle number, then this prescription is somewhat more complicated |19l . 

®The importance of surface emission has been stressed also in Refs. |23| for opaque sources which emit particles 
from a surface layer of finite thickness, then {Ar^^t}p < (Ar^jjjg)p reducing thereby Rout/ Rside ratio. 
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at early times happens at the periphery of the volume occupied by quark-gluon matter. 



2 Hydro-inspired parametrization of continuous hadronic freeze- 
out in central RHIC Au+Au collisions 

We assume that particles are emitted from some 3D hypersurface when the fireball is in a locally 
equilibrated (leq) state, 

/. = fie,A^,p) = [2.)-' (exp ( ^-^^^^^^;^^-^"^ ) ± (2) 

characterized by the temperature T{x) (which is common for all the particle species), collective 
velocity field Ufj_{x) and the chemical potentials ^i{x) depending on the particle species i. We 
refer to the part of 3D hypersurface with a time- like normal vector as v- (volume) freeze-out, or 
final breakup, and the rest of the hypersurface with a space-like normal as s- (surface) freeze-out, 
corresponding to continuous emission prior to the final breakup. We assume that particle densities 
are uniform at each v— and s— hypersurfaces separately, but, in general, are different from each 
other. Then the temperature and chemical potential in Eq. ^ are replaced with T^, Tg and ^u^^j, 
IJLs,i for the V— and s— freeze-out sectors, respectively. For the sake of simplicity and reduction of 
the number of parameters we do not interpolate smoothly the values of temperature and chemical 
potential between the volume and surface sectors of the hypersurface. A longitudinally boost- 
invariant expansion is a rather good approximation in mid-rapidity region at RHIC j24j, therefore 
we assume boost invariance for longitudinal hydrodynamic velocities in this region, vl = z/t, which 
allows us to parameterize t and z at 3D hypersurface a as ty^s = '''v,s{t^) cosh(r/), z^^s = 'Tv,s{t^) sinh(7/) 
in both sectors of the freeze-out hypersurface, where r = \Jt'^ — z'^ is the Bjorken proper time j25j 
and T] is the longitudinal fluid rapidity, r\ = tanh~^f2,; in addition, r is the absolute value of 
transverse coordinate rx = (3;,y) = (r cos r sin 0). Taking into account the transverse velocity 
component f (r, r) in the longitudinally co-moving frame, we obtain for the collective 4-velocity 

g (r, r/) = (cosh ry cosh r/!^'** , sinh ry^'* cos 0, sinh r/!^'** sin 0, sinh r/ cosh rjl^^), (3) 

where rjl^^ {r,T^^s) is the transverse fluid rapidity, tjt = tanh""*^ w(r, r). The particle 4-momentum 
can be expressed through the momentum rapidity y, transverse momentum px and transverse mass 

rriT = \J m? + p^, 

= {niT cosh y , Pt , sinh y ) . (4) 

Then ^ 

= cosh (y - r/) cosh r/T - sinhr/T. (5) 

The volume and surface elements of the freeze-out hypersurface cr(x) take the form 

da'"^^ = v^li^da = ±Ty^s{r)dr]drxdry{coshri, -dTy^s/droo, -dr^^s/dry, - sinhr/), (6) 

where 

V s _ (coshr/, -dr^^s/dr^, -dr^^^/dr^, - sinh??) 

~ (±(1 - {dT,,s/dr,r - {dT,^s/dryY)Yl^ ' ^ ' 
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and "+" and "-" correspond to the v— and s— sectors of the hypersurface respectively.'' The 
surface "moves" with a velocity less than the speed of light, |drs(r)/dr| < 1, hence \dTs{r) / dr\ > 1. 
If the 2D surface of the 3D s— freeze-out sector "moves" outwards, cLts/ drx^drg/ dry > 0, then the 
product p^n^ will be negative if the pT of particles are sufficiently small. In that case such particles 
cannot really escape from the system. For particles with fairly high transverse momenta one has 
p^Ufj, > 0, and they will be emitted from the surface of the system into vacuum and stream freely 
towards detectors. To take this into account we use a modified Cooper-Frye prescription with the 
substitution [HI HO] 

fl.eq.{x,p) e{p^n^'{x))fl,eqix,p) (8) 

to get the single particle spectra pod^N/d^p for diverse particle species and the correlation function 
of pions C{p, q): 

d^N f 

Po-^ = J P''d(Tf,9{p^n''{x))fieqXx,p) (9) 

C{p q) = l + I Jp^'da^e{p^n^'{x))fl,eq.{x,p)exp{iqx)\'^ ^^^^ 

{f Pida^e{p'^nf,{x))fieq.{x,pi)) {j P2da^e{p'^n^{x))fi,eq.{x,p2)) 

Here p = {pi + P2) l"^-, q = Pi — P2- To calculate the correlation function in the region of correlation 
peak we utilize the smoothness mass shell approximation: 

, , , \ ! P*^da^e{p*n^{x))fi,eq.{x,p*)&Mnx)? 
C(p, gjssH (11) 

Up*>'da^9{p*^^n^{x))fi,,q,{x,p*)y 

where p*^' = m? + ^^). 

We assume that the source of particles in central Au + Au RHIC collisions has cylindrical 
symmetry and is strongly stretched out in the beam direction {long— direction): the longitudinal 
size is much larger than the corresponding length of homogeneity [261 127j . The latter allows us to 
neglect the finiteness of the system in the longitudinal direction when we calculate particle spectra 
in the mid-rapidity region around y « 0. As for the transverse direction, we assume that the source 
has a finite geometrical size encoded in the limits of integration over r: < r < Rf for t>— emission 
and Ri < r < Rf for s— emission, where Ri and Rf are the initial and finite effective radii of the 
system. 

To fix a model one needs to define T(r) and the transverse rapidity r]T{r). We choose the simplest 
ansatz, aimed at catching the main features of particle emission. We suppose that the volume 
(space-like) emission happens as assumed in the blast-wave model at the constant Ty{r) = Tf, 
and the surface emission takes place during the whole time of hydrodynamic evolution of system: 
Tj < Ts < Tf, where Tj is the initial time of thermalization. Because the fireballs formed in RHIC 
collisions expand in transverse direction with rather high transverse velocity, we can assume that 
dTs{r)/dr > for the surface freeze-out.* For the sake of simplicity we parameterize Ts{r) by a 
linear function, 

Ts{r) = a-r + b. (12) 

^Note that the above representation of elements of the freeze-out hypersurface a{x) in the s-sector is, of course, 
not unique and is used because of its convenience for the chosen type of (j{x) (see below). It should be changed if 
drs{T)/dT = on some part of this hypersurface. 

^This is in full correspondence with the form of the hypersurface of constant densities (particle number and energy) 
in analytic ellipsoidal solutions of relativistic hydrodynamics I28| . It seems that initial or fast developed transverse 
velocities are required in order to provide a positive t — r correlation in s- sector of freeze-out hypersurface. 
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Taking into account that Ts{Ri) = Ti and Ts{Rf) = Tf, where Ri and Rf are the initial and final 
transverse radii respectively, Ri < Rf, we obtain that 

- (13) 



Rf-Ri 
Tf-n 



Rf ~ Ri 



R^. (14) 



Please note that the chosen form of Ts{r) is just an approximation that reflects the main properties 
of the surface freeze-out. For example, the chosen constant velocity of outward motion of the 
transverse surface, 1/a, is associated with the mean (averaged over time) velocity of the surface, 
which generally can move with time-dependent velocity. 

The form of freeze-out r(r) is presented in Fig. 1. It is convenient to parameterize the volume 
as well as surface sectors of the 3D hypersurface by tt and longitudinal rapidity ??, therefore the 
integration in Eqs. ©, (|TT|) is carried out separately for those sectors: / d(T^[...] J da^[..]y + 
fda;[..]s. 

To specify the transverse rapidity r]T{r), we first connect smoothly the values of transverse 
rapidity in the volume and surface sectors of the hypersurface on their "border". Second, since a 
linear r dependence of the radial velocity profile was found in hydrodynamic simulations j2S] for 
small transverse collective velocities v ^ rjT near the center of the system at fixed r, we assume 
that the transverse flow rapidity profile depends linearly on r for volume freeze-out at t = Tf, 

vUr,rf)=v¥^'^, (15) 
Rf 

where rjjp"'^ is the maximal transverse rapidity of the fluid. As for the transverse rapidity profile in 
the surface sector of the freeze-out hypersurface, we analyzed various parameterizations. Fits to 
the data prefer a non-zero initial transverse velocity at Tj = 1 fm/c, which is physically reasonable. 
Indeed, it seems to be natural that the transverse flow velocity at freeze-out surface at any time 
(or, at least, on average) should not be smaller than the average "velocity" of the surface, 1/a. 
Such an "initial" velocity, VT{Ti = 1 fm/c), can be developed even at the non-thermal partonic 
stage before Tj. Also a non-zero initial velocity does not contradict an approximate linearity with 
radius of rapidity profile at early times, in particular at Tj. One of the simplest parameterizations 
which satisfies these properties and allows us to connect smoothly the values of transverse rapidity 
in the volume and surface sectors is the following transverse rapidity distribution: 



(Tf - + R) 

Note that Ts{Rf) = Tf and r]?^{Rf ,T,{Rf)) = r]^{Rf,Tf) = r/?^'^^. 

In the next Section we apply our model to description of pion, kaon and proton transverse 
momentum spectra and pion interferometry radii in mid-rapidity in RHIC Au+Au central collisions 
at ^/sj^ = 200 GeV. 

3 Results and discussion 

We assume a common enclosed freeze-out hypersurface for all particles species. Then the model 
contains 7 parameters: Tj, Tf, Ri, Rf, Ty, Tg, rjjp"'^. As commonly accepted, we assume that at 
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Ti = 1 fm/c the system is already in local thermal equilibrium, and so we will not consider this 
value free parameter. 

The values of chemical potentials /^sj, fJ^vj which are different for different particle species j, 
are not free parameters and are responsible for the absolute normalization of the spectrum. In 
details, the situation is somewhat more complicated since resonance decays contribute to both the 
particle spectra and interferometry radii j30|. About 1/2 — 2/3 of pions, kaons and protons come 
from those decays after chemical freeze-out happens in high-energy heavy ion collisions j^H 1321 127j . 
The corresponding contribution to pion spectra and interferometry radii was estimated in Ref. j33j 
within the hydrodynamic model developed in |22j . According to the results of Ref. [HHI , resonance 
decays do not change significantly the slopes of the transverse pion spectra in the region of interest 
for kinetic freeze-out temperatures T„ and flows at RHIC collisions. Note that the resonance 
emission from the surface is negligible since heavy resonances have too small heat velocities to 
escape from the expanding system. As for the pion interferometry radii, the decays of short- 
lived resonances result in an augmentation of them as compared with the picture of pure thermal 
emission, especially in the region of small pT- This resonance contribution to the pion interferometry 
radii decreases with the increase of transverse momenta of the pairs and, as it was found in Ref. 
PH] . is smaller than 6 — 14 % for pT > 0.3 GeV (see also 1311 )• 

Based on these studies, we neglect the changes of slopes of particle transverse spectra due to 
resonance decays and fit the pion interferometry radii starting from niT ~ 0.3 GeV where the 
resonance contributions are fairly small. Finally, to fix fJ-vj, fJ-sj for thermal (primary) particles at 
some Tjj and Tg, one needs to know what fraction of measured particles is produced by resonances 
after thermal (kinetic) freeze-out. Such a problem is not fully solved until now and needs further 
analysis. Since the main aim of this paper is to describe the main features of the physical picture 
of particle emission in central RHIC Au+Au collisions rather than to fine tune the parameters and 
perfectly describe the data, we fit the fraction of thermal (direct) particles within reasonable limits 
discussed above; our final choice is presented in the last column of Table 1. 

A preliminary analysis shows that a good fit to spectra and correlations formed on an enclosed 
hypersurface cannot be obtained if the surface temperature is taken within the standard range 
for kinetic freeze-out: 100 MeV < < ttItt; it should be taken higher than the pion mass m,r-^ It 
forces us to conclude that the temperature of surface emission is close to the chemical freeze-out 
temperature, or the hadronization one if chemical freeze-out happens at about the same time, as 
it is argued in Ref. [SHI- Since the periphery of an expanding system is more dilute, the matter 
there is in the hadronic phase and hadrons can easily escape from this layer into the surrounding 
vacuum despite the relatively high temperature; this was argued in detail in Ref. ^Hl- This natural 
picture of surface emission at earlier times < r/ gives us a possibility to assume i^s = for pions 
and kaons, accounting for chemical equilibrium of pions at chemical freeze-out (see, e.g., |321) and 
almost zero net strangeness in central RHIC Au+Au collisions at the top energy. The results of 
our fit, which are presented in Tables 1 and 2, demonstrate that the best fit fixes the value of Tg at 
150 MeV, which is indeed very close to the chemical freeze-out temperature 157 it 3 MeV obtained 
in Ref. j37j . This justifies our assumption a posteriori. The best fit to proton transverse spectra 
is achieved at //s^p = 40 MeV, which is slightly more than the value fiB = 28.2 it 3.6 MeV found 
in Ref. [HZI for the baryochemical potential. ^'^ Concerning the values /U^j, one readily sees that 

®The same conclusion that the best fit of the pion HBT radii is obtained when the temperature of surface emission 
from a hypersurface with a spacehke normal is significantly higher than the temperature of volume emission is made 
also in Ref. [35|, where SPS data for 158 AGeV central Pb+Pb collisions are analyzed. 

^"Note that it is not only the absolute values but also the spectrum slopes that are affected by the values of f^sj 
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they are not free parameters of the model, since they should be fixed by the measured respective 
rapidity densities, dNj/dy. 

The best fit to data is obtained for the parameters shown in Tables 1 and 2. We also present 
there parameters that were kept constant through the fit procedure. As one can see from Figs. 2 — 4 
and Tables 1, 2, that good description of pion, proton and kaon spectra and pion interferometry 
radii is achieved in the mid-rapidity region of central -^/s/vTV = 200 GeV Au+Au collisions. The 
STAR collaboration data are presented for illustration. We did not fit these data because they are 
limited to relatively small pT values as compared to the PHENIX data, whereas the surface emission 
gives a significant contribution in the region of high transverse momenta. It is noteworthy, that the 
best fit is achieved when the system already at = 1 fm/c expands in the transverse direction with 
a collective velocity field that reaches at the system border (at r = Ri) the value v{Ri,Ti) 0.34 
c. It may indicate that a noticeable contribution into buildup of transverse flow at RHIC energies 
is given at times earlier then 1 fm/c, possibly pre-thermal stages of evolution. It is interesting 
to note that hydrodynamic simulations of heavy ion collisions at RHIC energies also indicate 
the necessity of an initial (pre-hydrodynamic) transverse flow to better account for slopes of the 
observed spectra. 

TABLE 1. 

Model parameters from fits of Au+Au ^/snn — 200 GeV data on tt^ , K~ and p single-particle spectra and 
TT^ HBT radii. The ratio of detected particle numbers to thermal (primary) ones (last column) as well as 
fis for pions and kaons were kept constant during fit procedures. 





fly MeV 


/xs MeV 


(iAf=9 /dy 
dm'^/dy 


TT 


53 





2.2 




45 





2.2 


p 


280 


40 


3.5 



TABLE 2. 

Model parameters from fits of Au-|-Au ^/snn — 200 GeV data on tt", K~ and p single-particle spectra and 
TT^ HBT radii. The value of was kept constant through fit procedures. 





Ti fm/c 


Tf fm/c 


Ri fm 


Rf fm 


Ty MeV 


Ts MeV 


'IT 


TT , K , p 


1 


10 


7 


11.3 


110 


150 


0.73 



The solid lines in Figs. 2 — 4 demonstrate the fit according to our model, and the dashed 
lines indicate the contributions to particle spectra and interferometry radii from volume freeze-out 
solely. The latter corresponds, in fact, to blast-wave model calculations with our volume freeze- 
out parameters and Ar = 0. The experimental data are taken from Refs. |39| I40| I41j . We 
used Boltzmann approximation for the distribution functions in the single-particle spectra and 
interferometry radii calculations for all particle species. One can see from Fig. 2 that surface 
emission leads to a noticeable enhancement of the single-particle momentum spectra in high px 
region, resulting thereby in an increase of the inverse spectra slope, or effective temperature. This 
happens because the particles with relatively low px are not emitted from the surface part of the 
freeze-out hypersurface that moves outwards and so their contribution to the spectra is suppressed, 

because chemical potentials determine what part of particles are emitted from the surface, and these particles have, 
of course, different spectrum slopes than the ones emitted from the volume part of the freeze-out hypersurface. 
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while the particles with sufficiently high px are liberated freely. This effect is taken into account, 
as discussed in the previous Section, by means of multiplying the distribution functions by the 
step function, see Eq. Q, thus restricting the freeze-out to those particles for which the product 
p^n^ is positive. Note that in the blast-wave parametrization, which takes into account only volume 
emission, such a high effective temperature can only be obtained at the cost of a very high transverse 
flow {r]jp°-^ 0.9 and so v^"-^ ^ 0.7 c, while in our parametrization r/™"^ 0.72, v'"^"'^ 0.6 c) at 
freeze-out |22j . which seems less realistic and more difficult to reproduce in hydrodynamic models 
of the system evolution. 

The Figures 3, 4 show our results for the interferometry radii of negative pions. Note that the 
interferometry radii are measured by the PHENIX Collaboration for — 30 % centrality events 
at ^/sNN = 200 GeV 41 . Therefore, to adjust the results of the PHENIX Cohaboration to the 
most central — 5 % centrality bin, we increase the corresponding interferometry radii by a factor 
calculated in accordance with the Npart dependence of the Bertsch-Pratt radius parameters found 
in Ref. [l^. As one can see from Fig. 3, the interplay between (i) a long emission time, r/ — r,, 
which gives an increase of the Rout radius, and (ii) positive ro^t — t correlations at the surface part 
of the hypersurface where drg/dr > 0, which reduces the radius, provides a moderate increase of 
Rout as compared with the volume emission at constant tM The Rgide radius is almost not affected 
by the surface emission. As to Riong: it is less than in the volume emission case. The reason 
is that Riong is strongly affected by the value of the emission time r j^Hl and, thereby, is reduced 
by the emission from early times t < tj. This can explain an apparently low freeze-out value of 
T in relativistic heavy ion collisions if it is extracted assuming instantaneous volume emission at 
constant r. 

Figure 4 represents our result for the Rout/ Rside ratio. One can see that the volume blast-wave- 
like emission (at a constant r) leads to a rather low Rout/ Rside ratio. The "extra" reduction of 
Rout is typical for a sharp freeze-out hypersurface with step density profile in the out-direction, 
n{r) ~ e{R 

~ Tout) and the maximal velocity on the edge. The same boundary effect of the 
additional broadening of the correlation function at large momenta in the case of sharp step-like 
freeze-out in the long-direction have been found and described in detail in Ref. 42_ . As one can see 
from Fig. 4, this effect becomes softer for an enclosed freeze-out hypersurface despite the the fact 
that maximum velocity is also achieved at maximum transverse radius which the system occupies 
during its evolution. However, negative t — rout correlations can totally destroy the effect of relative 
reduction of Rout- That is why the positive t — rout correlations in our model are very important 
to achieve a satisfactory description of the experimental data. 

4 Summary 

We have studied the effect of continuous particle emission on particle spectra and interferometry 
radii, based on hydro-inspired parametrization of the enclosed freeze-out hypersurface and distri- 
bution functions. Our model generalizes the "blast-wave" parametrization, which has difficulties 
in the fitting of the HBT radii and leads to too large transverse velocities. We demonstrate that 
7r~, K~ and p single particle spectra and vr" interferometry radii at RHIC can be described by 
a physically reasonable set of parameters when a rather protracted emission (~ 9 fm/c) from the 
surface sector of the hypersurface that moves "outward" is taken into account. The phase-space 

^^The conclusion that a long-lived fireball can be compatible with HBT measurements is made also in Ref. |35| 
^^An independence of Rside on transport opacity was also found in a covariant transport model |18|. 
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distribution of particles emitted from the surface prior to the fireball breakup is similar to the 
local chemical equilibrium distribution with parameters rather close to the ones found at chemical 
freeze-out in RHIC Au+Au collisions. This could indicate a transformation of quark-gluon matter 
to hadron gas at the periphery of the system with subsequent hadronic emission just after such a 
transition from the boundary of the fireball during the whole time of the system's hydro evolution. 

Thereby, the simple physical arguments requiring the freeze-out hypersurface to be enclosed for 
finite expanding systems result in a successful description of hadronic spectra and interferometry 
radii. This means that particle emission from the early stages of fireball evolution has to be taken 
into account in advanced hydrodynamic models of A+A collisions. A quantitative description of 
a continuous escaping of particles can be made, in particular, within the hydro-kinetic approach to 
spectra formation proposed in Ref. 
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Figure 1: The enclosed freeze-out hypersurface T(r). See the text for details. 
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Figure 2: Comparison of the single-particle momentum spectra measured by the PHENIX Col- 
laboration with the model calculations performed for whole enclosed hypersurface (solid lines) and 
for shelf-like part, see Fig. 1, solely (dashed lines). For convenience of plots location the measured 
p spectra are reduced in 100 times and K~ spectra in 8 times. The STAR data are not fitted and 
are presented for illustration. 
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Figure 3: Comparison of the pion source Rout, Rside, Rlong radii measured by the STAR (stars) and 
PHENIX (boxes) Collaborations with model calculations performed for whole enclosed hypersurface 
(solid lines) and for shelf-like part only (dashed lines). The PHENIX data are adjusted for most 
central bin (see the text for details). 
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Vs" = 2 00 GeV; tt"; * STAR; ■ PHENIX 
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Figure 4: Comparison of the measured by the STAR (stars) and PHENIX (boxes) Collaborations 
pion Rout/Rside ratios with model calculations performed for whole enclosed hypersurface (solid 
line) and for shelf- like part only (dashed line). 
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